Vehicular traffic flow at a non-signalised intersection 
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We have developed a modified Nagel-Schreckenberg cellular automata model for describing a 
conflicting vehicular traffic flow at the intersection of two streets. No traffic lights control the 
traffic flow. The approaching cars to the intersection yield to each other to avoid collision. Closed 
boundary condition is applied to the streets. Extensive Monte Carlo simulations is taken into 
account to find the model characteristics. In particular, we obtain the fundamental diagrams and 
show that the effect of interaction of two streets can be regarded as a dynamic impurity located at 
the intersection point. Our results suggest that yielding mechanism gives rise to a high total flow 
throughout the intersection especially in the low density regime. In some ranges of densities, yielding 
mechanism even improves and regulates the flow in comparison to the absence of perpendicular flow. 
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I. INTRODUCTION 



Modeling the dynamics of vehicular traffic flow has 
constituted the subject of intensive research by statis- 
tical physics and applied mathematics communities dur- 
ing the past years |E IE IE IE Q- In particular, cellu- 
lar automata approach have provided the possibility to 
study various aspects of these truly non-equilibrium sys- 
tems which still are of current interest H|- Besides 
various theoretical efforts aiming to understand the ba- 
sic principles governing the spatial-temporal structure of 
traffic flow, considerable attempts have been made to- 
wards realistic problems involving optimization of vehic- 
ular traffic flow. While the existing results in the con- 
text of highway traffic seem to need further manipula- 
tions in order to find direct applications, researches on 
cit y tr affic have more feasibility in practical applications 
[E EE El EE EE EI EE EE El ■ We believe that opti- 
misation of traffic flow at a single intersection is a sub- 
stantial ingredient for the task of global optimisation of 
city networks [IE]. Isolated intersections are fundamen- 
tal operating units of complex city networks and their 
thorough analysis would be inevitably advantageous not 
only for optimisation of city networks but also for local 
optimization purposes. Recently, physicists have paid 
notable attention to controlling traffic flow at intersec- 
tions and other traffic designations such as roundabouts 
EE HE HE EE EE EI EE HE HE Hi- In this respect, 
our objective in this paper is to study another aspect 
of conflicting traffic flow at intersections. In principle, 
the vehicular flow at the intersection of two roads can 
be controlled via two distinctive schemes. In the first 



scheme, which is appropriate when the density of cars in 
both roads are low, the traffic is controlled without traf- 
fic lights. In the second scheme, signalized traffic lights 
control the flow. In the former scheme, approaching car 
to the intersection yield to traffic at the perpendicular di- 
rection by adjusting its velocity to a safe value to avoid 
collision. According to driving rules, the priority is given 
to the nearest car to the intersection. It is evident that 
this scheme is efficient if the density of cars is low. When 
the density of cars increases, this method fails to opti- 
mally control the traffic and long queues may form which 
gives rise to long delays. At this stage the implementa- 
tion of the second scheme i.e.; utilizing traffic lights is 
unavoidable. Therefore it is a natural and important 
question to find out under what circumstances the in- 
tersection should be controlled by traffic lights? More 
concisely, what is the critical density beyond which the 
non-signalised schemes begins to fail. In order to capture 
the basic features of this problem, we have constructed a 
cellular automata model describing the above dynamics. 
This paper has the following layout. In section II, the 
model is introduced and driving rules are explained. In 
section III, the results of the Monte Carlo simulations are 
exhibited. Concluding remarks and discussions ends the 
paper in section IV. 



II. 



DESCRIPTION OF THE PROBLEM 
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We now present our CA model. Consider two perpen- 
dicular one dimensional closed chains each having L sites. 
The chains represent urban roads accommodating unidi- 
rectional vehicular traffic flow. They cross each other at 
the sites i\ = %2 = ■§ on the first and the second chain 
respectively. With no loss of generality we take the direc- 
tion of traffic flow in the first chain from south to north 
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FIG. 1: Intersection of two uni-directional streets. 



and in the second chain from east to west (see Fig.l for 
illustration) . The discretisation of space is such that each 
car occupies an integer number of cells denoted by L car . 
The car position is denoted by the location of its head 
cell. Time elapses in discrete steps of At sec and veloc- 
ities take discrete values 0, 1, 2, • • • , v max in which v max 
is the maximum velocity of cars. 

To be more specific, at each step of time, the sys- 
tem is characterized by the position and velocity con- 
figurations of cars. The system evolves under the Nagel- 
Schreckenberg (NS) dynamics (29|. Let us briefly explain 
the NS updating rules which synchronously evolve the 
system state from time t to t + 1. We denote position, 
velocity and space gap of a typical car at timestep t by 



(0 



«''' and gW respectively. The same quantities for 
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its leading car are correspondingly denoted by x 

and gf . We recall that gap is defined as the distance 
between front bumper of the follower to the rear bumper 
of its leading. More precisely, g(t) = xi(t) — x(t) — L car . 
Concerning the above considerations, the following up- 
dating sub steps evolve the position 
and the velocity of each car in parallel. 

1) Acceleration: 

v (t+l/3) ._ min (y(t) + l )tw ) 

2) Velocity adjustment : 

v (t+2/3) ._ min ( 5 (t+l/3) )V (t+l/3)) 

3) Random breaking with probability p: 

if random < p then v^ t+1 ^ :— max(v^ t+2 ^' i ' 1 — 1, 0) 

4) Movement : := x^ + v (t+1) 

The yielding dynamics in the vicinity of the intersec- 
tion is implemented by introducing a safety distance D s . 
The approaching cars (nearest cars to the crossing point 



FIG. 2: Two approaching cars to the intersection yield to 
each other provided their distances to the crossing point are 
both less than a safety distance. 



i = U) should yield to each other if their distances to 
the crossing point, denoted by d\ and di for the first and 
second street respectively, are both less than the safety 
distance D s . In this case, the movement priority is given 
to the car which is closer to the crossing point. This car 
adjust its velocity as usual with its leading car. On the 
contrary, the further car, which is the one that should 
yield, brakes irrespective of its direct gap. The simplest 
way to take into this cautionary braking is to adjust the 
gap with the crossing point itself. This implies that the 
yielding car sees the crossing point as a hindrance. In this 
way, the model is collision-free. Figure two illustrates the 
situation. 

Let us now specify the physical values of our time and 
space units. Ignoring the possibility of existence of long 
vehicles such as buses, trucks etc, the length of each cell 
is taken to be as 4.5 metres which is the typical bumper- 
to-bumper distance of cars in a waiting queue. There- 
fore the spatial grid Ax equals to 4^- m. We take the 
time step At = 1 s. Furthermore, we adopt a speed- 
limit of 75 km/h. The value of v max is so chosen to give 
the speed-limit value 75 km/h or equivalently 21 m/s 
. In this regard, v max is given by the integer party of 
21 x L car /4.5. For instance, for L car = 5, v max equals 
23. In addition, each discrete increments of velocity sig- 
nifies a value of Av = L 4 5 m/s which is also equivalent 
to the acceleration. For X car = 5 this gives the comfort 
acceleration 0.9 Moreover, we take the value of ran- 
dom breaking parameter at p = 0.1. In the next section, 
the simulation results of the above-described dynamics is 
presented. 



III. 



MONTE CARLO SIMULATION 



The streets sizes are equally taken as L\ = = 
1350 m and the system is update for 10 6 time steps. 
After transients, two streets maintain steady-state cur- 
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FIG. 3: Ji vs pi for various values of p2. 
icar = 5 cells. The road length is 1350 metre. 



25 m and FIG. 4: J2 vs pi for various values of p2- Ds = 25 m and 
icar = 5 cells. The road length is 1350 metre. 



rents, defined as the number of vehicles passing from a 
fixed location per a definite time interval, denoted by 
Ji and J2 . They are functions of the global densities 
pi _ Ni x L ca .r and p2 _ Jv 2 x£ cai . w h ere Nl and N 2 are 

the number of vehicles in the first and the second street 
respectively. We kept the global density at a fixed value 
P2 in the second street and varied pi . Figure (3) exhibits 
the fundamental diagram of the first street i.e.; J\ versus 
Pi- 

It is observed that for small densities P2 up to 0.05, 
Ji rises to its maximum value, then it undergoes a 
short rapid decrease after which a lengthy plateau re- 
gion, where the current is independent of p\, is formed. 
Intersection of two chains makes the intersection point 
appear as a site-wise dynamical defective site. It is a 
well-known fact that a local defect can affect the low 
dimensional non-equilibrium systems on a global scale 
[I], ED, [II [H HI M, H M, [Si. This has been con- 
firmed not only for simple exclusion process but also 
for cellular automata models describing vehicular traf- 
fic flow [39|, Hfj. Analogous to static defects, in our 
case of dynamical impurity, we observe that the effect 
of the site-wise dynamic defect is to form a plateau re- 
gion p G p + ] in which A = p + — p_ is the extension 
of the plateau region in the fundamental diagram. The 
larger the density in the perpendicular chain is, the more 
strong is the dynamic defect. For higher p2, the plateau 
region is wider and the current value is more reduced. 
After the plateau, J\ exhibits linear decrease versus pi 
in the same manner as in the fundamental diagram of a 
single road. In this region which corresponds to p\ > p+ 
the intersecting road imposes no particular effect on the 
first road. Increasing P2 beyond 0.05 gives rise to sub- 
stantial changes in the fundamental diagram. Contrary 
to the case P2 < 0.05, the abrupt drop of current after 
reaching its maximum disappears for P2 > 0.05 and J\ 
reaches its plateau value without showing any decrease. 
The length and height of the plateau does not show sig- 
nificant dependence for P2 G [0.05,0.8]. This marks the 



efficiency of the non signalised controlling mechanism in 
which the current of each street is highly robust over the 
density variation in the perpendicular street. When p 2 
exceeds 0.8, the plateau undergoes changes. Its length 
increases whereas its height decreases. We now consider 
the flow characteristics in the second street. Although 
the global density is constant in street 2 its current J2 is 
affected by density variations in the first street. In Figure 
(4) we depict the behaviour of J2 versus pi. 

For each value of p2, the current J 2 as a function of 
pi exhibits three regimes. In the first regime in which 
pi is small, J2 is a decreasing function of p±. After- 
wards, J2 reaches a plateau region (second regime) which 
is approximately extended over the region pi G [0.1, 0.8]. 
Eventually in the third regime, J2 exhibits decreasing 
behaviour towards zero. Analogous to Ji, the existence 
of wide plateau regions indicates that street 2 can main- 
tain a constant flow capacity for a wide range of density 
variations in the first street. The other feature is that in 
fixed pi, J2 is an increasing function for small values of 
P2- This is natural since the current in street 2 has not 
reached its maximal value. This increment persists up to 
P2 = 0.05. Beyond that, for each pi, J 2 saturates. In the 
plateau region, the saturation value is slightly above 0.16. 
The current saturation continues up to P2 — 0.8 above 
which J2 again starts to decrease. We note that the be- 
haviours depicted in Ji — pi and J2 — pi diagrams are 
consistent to each other. Due to the existence of 1 ^ 2 
symmetry, the J2 — P2 diagram is identical to J\ — pi and 
J\ — P2 is identical to J2 — pi- In order to find a deeper 
insight, it would be illustrative to look at the behaviour 
of total current J to t = J\ + J2 as a function of density in 
one of the streets. Fig. (5) sketches this behaviour. 

For p2 < 0.05, the maximum of Jtot lies at pi = 0.1. 
However, for P2 > 0.05, the maximum shifts backward 
to pi = 0. According to the above graphs, after a short 
increasing behaviour, J tot enters into a lengthy plateau 
region. Evidently for optimisation of traffic one should 
maximize the total current Jtot- The existence of a wide 



4 



"T" 



"T" 



A ^ & H 
.'; AAA. 



A A-'X. ; 



□ (I, = 0.01 

A Pj = 0.02 

V % = 0.03 

Is p 2 =0.04 

a—- p= 0.05-0 

\\ - -0- - P 2 = 0.S5 

V V V V V V V V V i; S - -O (^ = 0-9 

''V > >-> >■ > > > > 

^..^.^_ < ]_. < _..^.. <] ..^ ] ._ C 3_ <1 _.^_. 



- -© - O -O -G-O- o -© ■ e- -o-o-o- e- -e- -O -Q-Ov 




J 1 02 

0.15 



4> o 



-n— n- b- a- b -e -b -b -b-b- a. 

□ Lcar = 1 \ 

O Lear = 3 \ 

O Lear = 5 R 
A Lear = 7 \ 

-O O O O O G-O G O O-Q-O-Q-h 

<> <> <-> o o o o o o O O 



■A A A -A— A— A— A A A A -A A A-A-AA 

R, = 0.5 



\ 



0.25 



0.5 



0.75 



FIG. 5: Total current J to t vs pi for various values of p2. 
D 3 — 25 m and L car = 5 cells. 



FIG. 7: Ji vs pi for various L car with p2 = 0.5. 
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FIG. 6: Ji vs pi for various L car . Global density of the 
second street is kept fixed at p2 = 0.05. 

plateau region in J to t suggests that yielding mechanism 
can be regarded as an efficient method in the plateau 
range of density in the first street. Let us now consider 
the role of L car . The cellular nature of our model permits 
us to adjust the cell length Ax in such a way to reproduce 
a reasonable acceleration. Our simulations demonstrate 
that currents exhibit significant dependence on Ax or 
equivalcntly on L car . This is exhibited in figures (6,7). 

While the structure of the fundamental diagram does 
not qualitatively change, the values of J\ notably depend 
on L car . For both p 2 = 0.05 and 0.5, J\ is a decreasing 
function of L car . The reason is that larger L car gives rise 
to higher acceleration. 

Analogous to J\, the dependence of J2 on L car is con- 
siderable as is shown in figures (8,9). Variation of L car 
does not lead to change the generic behaviour but rather 
changes the current values. Due to the same reason which 
was explained, smaller L car gives higher currents. In the 
case L car — 1, the transition of Ji from the plateau re- 
gion to the linear decreasing segment is much smoother 
compared to the other values of L car greater than one. 
Since the currents in the plateau region do not depend 



FIG. 8: J2 versus pi for various L car . p2 = 0.05. 

on density, therefore the higher acceleration gives rise to 
larger currents. 

Finally, we have also examined the effect of varying 
the safety distance D s . Our simulations do not show any 
significant dependence on D s . This is due to unrealistic 
decceleration in the NS model. 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We have investigated the flow characteristics of a 
non signalised intersection by developing a Nagel- 
Schreckenberg cellular automata model. In particular, we 
have obtained the fundamental diagrams in both streets. 
Our findings show yielding of cars upon reaching the in- 
tersection gives rise to formation of plateau regions in 
the fundamental diagrams. This is reminiscent of the 
conventional role of a single impurity in the one dimen- 
sional out of equilibrium systems. The performance of 
non-signalised controlling mechanism is especially effi- 
cient when the car density is considerably low in both 
streets. The existence of wide plateau region in the total 
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scheme to the density fluctuations and offers an optimal 
method for controlling the traffic at low densities. Our 
CA model allows for varying space and time grids. By 
their appropriate adjusting, we are able to reproduce a 
realistic acceleration. In low densities, the currents ex- 
hibits notable dependence on the values of spatial dis- 
cretisation grid. Finally we remark that our approach is 
open to serious challenges. The crucial point is to model 
the yielding braking as realistic as possible. Empirical 
data are certainly required for this purpose. We expect 
the system characteristics undergo substantial changes if 
realistic yielding declaration is taken into account. 



FIG. 9: J2 versus p\ for various L car 
fixed in the second street at p2 = 0.5 . 



The density is kept 
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